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We find the field configurations and the propagation constants of the 
guided modes in a dielectric waveguide with curved axis and rectangular 
cross-section. Outside the guide, the refractive index is uniform. Inside, 
the index profile in the radial direction (intersection of the meridional plane 
and the plane of curvature) follows a parabolic law with the maximum at 
the center of the guide; in the direction perpendicular to the plane of curva- 
ture the index is either uniform or parabolic, again with the maximum at 
the center of the guide. The guide with mixed profiles has been proposed as 
an easy-to-support, low-loss, ribbon-like guide for millimeter and optical 
waves while the other, with parabolic profile in both directions, is similar to 
the "SELFOC®" or "GRIN" image transmitting guides. 

The axial field components are small compared to the transverse com- 
ponents and consequently the modes are almost of the TEM kind. Within 
the guide the field distribution along a quadratic profile is a parabolic 
cylinder function of order close to an integer, and is sinusoidal along the 
uniform profile. The field components outside of the guide decay almost 
with exponential law. 

Inside the SELFOC-like guide, the field distribution of the funda- 
mental mode is gaussian and except for the attenuation the characteristics 
of the beam are similar to those obtained for a guide in which the parabolic 
index profile is not truncated. 

The attenuation constant a of any mode is very sensitive to the radius of 
curvature R. Doubling R reduces a by several orders of magnitude. 

Fixing R and the difference of refractive index between the center of the 
guide and the edge of it, the attenuation constant a passes through a mini- 
mum for a guide width measured in the plane of curvature which is only a 
few beam-widths. 

Radiation loss for the fundamental gaussian mode is negligibly small if 
the distance between the center of the beam and the edge of the guide is two 
or more half beam-widths. 

Guides with rectangular index profile in the plane of curvature have less 
radiation loss than similar guides with truncated parabolic profile. 
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I. INTRODUCTION 

A dielectric guide in which the refractive index decreases with para- 
bolic law away from its axis acts as a lens-like medium. 1,2 The trans- 
mission through it is known even if the axis is not straight 3 and if the 
parabolic decrease is different in two orthogonal directions 4 (astigmatic 
guide) . 

Though extremely useful in many respects the parabolic medium is 
not realizable since it has ever-decreasing refractive index away from 
the axis and this in turn produces an untenable physical result. Thus 
though we know that in any realizable dielectric guide with curved 
axis, radiation losses are inevitable, 5 the modes in the parabolic medium 
with curved axis can have no radiation loss since the refractive index 
tending towards infinity far away from the axis prevents it. 

A more realistic model is achieved by truncating the parabolic index 
distribution. We begin, in Section II, studying the two dimensional 
guide, Fig. la, in which the index profile, Fig. lb, varies as a truncated 
parabolic function along the x axis and is independent of y while out- 
side of the guide the index is uniform. 

Later, this guide is modified in such a way that along y, the index 
profile is either rectangular, Fig. 2a, or another truncated parabolic 
function, Fig. 2b. 

The first of these guides has the index distribution of the dielectric 
thin-film guide proposed in Ref. 6 as a low-loss, easy-to-support ribbon- 
like guide for millimeter and optical waves. It has also the configuration 
of a possible guide for integrated optics. 7 This guide, with curved axis 
has been analyzed in Ref. 8 ignoring radiation due to curvature. In 
Section II, both the phase and attenuation coefficients of the guided 
modes are evaluated and compared to those in a similar guide with 
rectangular index profiles along both x and y. 

The results obtained for the guide with truncated parabolic profiles 
along x and y, Fig. 2b, are applicable, at least in order of magnitude, 
to "SELFOC" 9 or "GRIN" 10 fibers, and tubular gas lenses 11 with 
curved axes. 

Finally conclusions are drawn in Section III, while all the mathe- 
matics are given in the Appendix. 

II. MODES IN THE CURVED GUIDE 

Consider the two-dimensional curved guide in Fig. la. The parabolic 
refractive index within the guide is independent of y and equal to 

n,- = n[l - A(l + ^) 2 ] , (1) 
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Fig. 1 — (a) Two-dimensional truncated parabolic guide; (b) Refractive index 
profile; (c) Electric field distribution of the fundamental mode. 



where a is the width of the guide, n, the refractive index in the center 
of it and n(\ — A), the refractive index at the edges. Outside the guide, 
the index is 

n = n(l - A - A,). (2) 

We make the following assumptions: 

A « 1 



A, « 1 



and 



aV& 



«1 - 



4Aft 



(3) 



(4) 
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Fig. 2— (a) Inhomogeneous dielectric thin film guide; (b) "SELFOC®" or 
'GRIN" guides with rectangular cross-section. 



where A is the free-space wavelength and R the radius of curvature 
of the guide. The physical significance of inequality (3) is that the 
guided modes will have phase velocities quite comparable to that of a 
plane wave in a uniform medium of refractive index n. The inequality 
(4) insures that the amplitude of the field components at the edge of the 
guide are small compared to their maxima within the guide. In other 
words, most of the electromagnetic field is well confined within the 
guide, Fig. lc, and consequently the loss per wavelength is small com- 
pared to unity. Considering only guided modes with field configura- 
tions independent of y, we can group them in two families: TE and 
TM. The field components of any mode of the first family are E„ , 
H x and H s while those of the second are H„ , E x and E t . In each family 
the transverse components are far larger than the axial components 
and consequently both families are essentially of the TEM kind. 
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The transverse components E u , H x , H v and E x of both families 
have the same functional dependence within and without the guide. 
Therefore we will talk from now on of the E field meaning either one of 
those four components. 

Within the guide, and subject to the conditions (3) and (4), the E 
field distribution for the pth mode is essentially 



E = exp 



x + 2 " p 



w 



He, 



x + 2 - P 



w 



exp [i(k z z — <at)\ (5) 



in which the first two factors describe the field distribution along x, 
and the last gives the propagating wave dependence along the curvi- 
linear z axis. Similarly to the field distribution in the lens-like medium 
(a = oo ), the first factor is a gaussian with its maximum located at a 
distance 



P = 



8AR 



from the center of the guide. The normalizing 1/e half -width is 

w = J- 

* T 



a\ 



VsK 



(6) 



(7) 



The second factor in equation (5) is a Hermite polynomial of order p 
which is also centered at x = — (a/2) + p and the argument is normal- 
ized to w/2. Strictly speaking the expression (5) should have, instead 
of the Hermite polynomial, a Hermite function of order close to p. 
Interested readers can find the details in the Appendix. 

For the fundamental mode p = the Hermite polynomial is unity 
and the transverse field distribution is the well-known gaussian. 

The propagation constant k, = /3 + ia in equation (5) is complex 
and the phase and attenuation constants calculated in equations (36) 
and (37) are 

and 

2V^AdRp\ 



2 1 - M 

wkn K{\ + M). 



(8) 



exp^- 



01 



[«- 



d)> 



+*- 



p+i)(f 



2w 



- £s (i - <# 



(9) 
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in which 

i ~ 

0. = *»Vl - (±)'( P + i) , (10) 



wkn 



2 P 1 (a 



d = ±£ = ± « (11) 

a a (ft W ' 

(R = ^(2A)*fi, (12) 

A 

and the values of M and K can be found in equations (38) and (39). 
Let us discuss the physical meaning of some of these formulas. 

The phase constant given in equation (8) is the product of the 
phase constant ft^lO) of the lens-like medium with straight axis (R = 
a = oo ), multiplied by a bracket essentially equal to one; the two small 
terms contained therein take into account the curvature of the axis 
and the truncation of the parabolic profile. 

More interesting is the attenuation constant (9). The value a/2A Ra 
which is the normalized attenuation per radian has been plotted in 
Fig. 3 for the fundamental mode p = and A x = 0. The abscissa is 
the square of the guide width a normalized to the beam-width 2w or 
its equivalent (vna/X) Va/2" which is the guide width normalized to 
the free wavelength. The parameter used for the solid curves is the 
normalized radius of curvature (R(12). For a given radius of curvature 
the loss per radian is highly sensitive to the width of the guide and 
passes through a minimum at width 

a_ = fa\* 

2w \8/ ' 

For a wide range of values of (R, say 10 to 1000, that minimum loss 
occurs when the guide width is only a few beam-widths. 

The dotted lines are curves of constant d, that is constant ratio 
2 P /a between the beam displacement from the guide axis p and the 
guide half-width a/2. It is easy to understand the downward trend of 
these curves for large abscissas. Consider a guide with fixed geometry 
and decrease the wavelength X of operation. The beam remains at the 
same distance p from the guide axis but it becomes narrower and conse- 
quently the field at the edge of the guide and the radiation loss de- 
crease. It is surprising that the minimum radiation loss of the solid 
curves occurs when the beam displacement is a small part of the giduo 
width (d of the order of 0.1). 

Why do the solid lines have a minimum? For very narrow guides 
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Fig. 3 — Radiation loss in curved guides with truncated parabolic index profile. 
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(ct/2-w « 1), most of the electromagnetic field travels outside of the 
guide and any curvature of the axis introduces substantial radiation 
losses to this loosely guided beam. On the other hand, for very wide 
guides (a/2w >>> 1), any curvature of the axis displaces the beam close 
to one edge of the guide (cl close to unity) and once again substantial 
losses occur. There must be a minimum in between. 

It is interesting to compare the losses in these guides of truncated 
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parabolic index profile with guides of identical width but with rectangu- 
lar index profile of height nA. In Fig. 4, the solid curves are a repetition 
of some of those in Fig. 3, while the dotted ones have been reproduced 
from Ref . 12. The abscissa is again (a/2w) 2 which is identical to (ir/4)a/A 
in which 
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Fig, 4 — Radiation loss in curved guides with truncated parabolic index profile 
(solid curves) and with rectangular index profile (dashed curves). 
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is a dimension such that for a < A, the guide with rectangular index 
profile supports a single mode and for a > A, the guide is multimode. 

For the same radius of curvature, guide width, and same A on axis, 
the guide with truncated parabolic profile has more loss than the guide 
with rectangular profile. The difference is very marked for large abscissas, 
but this result should not be surprising because in the case of curved 
guides with truncated parabolic profile the beam travels close to one 
edge of the guide where there is little difference of refractive index 
between the inside and outside, while in the case of rectangular profile, 
though most of the power travels also close to one edge of the guide the 
full difference of refractive index nA is there to help in the guidance. 

In Fig. 5 we have plotted again the attenuation per radian as a 
function of (a/2w) 2 , but this time we use as parameter, the value of 



h = 



w 



which is the number of beam half-widths between the center of the 
beam and the external edge of the guide. The curves have asymptotes 
(dashed lines) parallel to both coordinates. 

For h ^ 2, A = 0.01, the attenuation per radian aR turns out to be 
smaller than 0.003, which is very small for most purposes. 

If the truncated parabolic profile is on a pedestal (A, 5^ 0), the 
losses are even smaller than those depicted in Fig. 4. The influence of 
A! in the attenuation constant (9), appears in the bracket of the ex- 
ponent. The other two terms are in general small compared to unity. 
Therefore even a modest value of A x , say Ai = A, is enough to reduce 
the losses depicted in Figs. 3 and 5 by several orders of magnitude. 

What happens when p ^ 0. From equation (9) we find as expected 
that for a given guide the radiation loss increases fast with the order p 
of the mode. The highest order mode that travels only slightly in- 
fluenced by the guide width is characterized by 



Pmax = 



w 



— A = h 2 — A 



Naturally p max is independent of A t , and when the beam center is close 
to a beam half-width from the edge, p max = 0. 

It is shown in the Appendix that if the refractive index profile along y, 
Fig. la, is not uniform but has either rectangular or truncated para- 
bolic shape, Figs. 2a and 2b, the guides have different phase constants 
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Fig. 5— Radiation loss in curved guides with truncated parabolic index profile. 
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than equation (8) but practically the same attenuation constant (9) 
provided that most of the electromagnetic field travels -within the 
guide. Therefore everything said about attenuation in this section 
applies to the three guides. 

For the following examples we will only use Figs. 3, 4 and 5 since all 
the important results and formulas are there. 
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2.1 Example A 

For a guide such that 

n = 1.5, 

A = 0.01, 

Ax = 0, 

a = 0.1 mm, 

X - 1 M , 

what is the radius of curvature R for which the loss per radian is of 
the order of 10~ 3 ? 

We calculate the abscissa and ordinate of Fig. 5 to be 

( a V _ rna /a _ 
\2w) " X \2 ~ A6 

and 

\/2AaK = 1.4- 10 -4 . 

The parameter h obtained from Fig. 5 is approximately 2 and we derive 

R = ^- = 3.9 mm. 
8Ap 

A very small radius indeed. 

2.2 Example B 

For integrated optics a guide with truncated-parabolic profile may 
have the following characteristics 

n = 1.5, 

A = 0.01, 

A, = 0, 

a = 10 M , 
X = 0.5m, 
R = 0.6 mm. 
What is the loss per radian? 
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From Fig. 3 or 4 we get the abscissa and parameter 

iS - «■ 

(R=^(2A)'^60. 
A 

Consequently the loss per radian results 

aR = 0.018. 

If instead of parabolic the index had been rectangular, from Fig. 4 
we deduce that the loss per radian would have been 0.00018, two 
orders of magnitude smaller. 

III. CONCLUSIONS 

For losses small enough, the field configurations and phase constants 
of the modes in dielectric guides, Figs. 2a and 2b, with curved axis 
and parabolic index profile on a pedestal, are quite comparable to 
those in a similar guide in which the parabolic profile is extended to 
infinity. 

The attenuation constant of a mode is very sensitive (exponential 
dependence) to the radius of curvature, size of the pedestal and order 
of the mode. The higher the order of the mode and the smaller the size 
of the pedestal the larger the loss. 

Quantitative results about the attenuation constant for the funda- 
mental gaussian mode in a guide without pedestal are given in Figs. 3, 
4 and 5 and in typical examples at the end at the preceding section. 
We find in these figures the loss per radian aR as a function of the guide 
width a, using as parameter the radius of curvature R, or the ratio 
between beam displacement p and guide width or the ratio between 
the beam distance from the edge of the guide, a/2 — p and the beam 
width w. The main conclusions are: 

(i) Doubling R reduces the attenuation constant a several orders 
of magnitude. 

(m) For any R, there is a guide width that minimizes the loss per 
radian. That dimension is only a few beam- widths. 

{Hi) For comparable characteristics, guides with rectangular pro- 
files have lower attenuation than those with truncated-para- 
bolic profile. Therefore if the transmission of images is not 
important, such as in the case of the ribbon-like guide of Ref . 6 
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and guides for integrated optics, rectangular index profiles are 
more attractive than parabolic profiles. 
(iv) The attenuation per 90° bend is smaller than 10" 3 in a guide 
such that the distance between beam center and the external 
edge of the guide is larger than a couple of half beam-widths, 
that is, if 

a _ 

2 



2w 



> 1 



APPENDIX 



Modes in Curved Guides 

With Truncated- Parabolic Index Profile 

We start studying the two-dimensional curved guide depicted in 
Fig. la in cylindrical coordinates. Later we will introduce a variation 
of the index profile along y. 

The parabolic refractive index distribution within the guide is 

r - & 



M-^)1 



n, ; =n 1 - A 1 +2 



(13) 



where a is the width of the guide, n the refractive index in the center 
and n(l — A) the refractive index at the edges. The refractive index 
outside the guide is 

n = n(l - A - A,)- (14) 

Assuming that the electromagnetic field does not vary along y and 
that the only component along that direction is H u , all the field com- 
ponents either inside or outside the guide are 

H.. = H 



E r m 



En = 



H 

we n,r 



dH 



coe nj dr 



r exp [i(v8 — oj<)] 



(15) 



where co is the angular frequency, e the refractive index of free space, 
and the indices i and o refer to the inside and outside of the guide. 
The resulting wave equation for both media is 

d*H 

dr 2 



r dr 



Vn\ -p)H = 



(16) 
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in which k = 2ir/X and X is the free space wavelength. Within the 
guide Tit is given by equation (13) and the wave equation can be re- 
duced to 

^ + b + h - ift + £o) 2 ]# = o (17) 

by making the following substitutions 



, _ 2(r - R) 
w 



&-|(l-4> f 



(18) 



in which 



, =» fc 2 K, (19) 

k 2 n t — k: 2 a 2 d I d\ 1 , Qft s 

'- A W -2^V-2j-2> (20) 



-a/?=a£S> (2i) 

, _ a 2 _ g = 2p ,__* 

fl i» 2 (R 4AiJ ' o ' k ; 

(R = ^(2A)*R, (23) 

A 



and 



A = £=■ (24) 

nV8A 

Furthermore, equation (17) has been derived making the following 
simplifying assumptions 

A«l, 

A 1 «l, (25) 

aVA 4 AS 

The physical significance of iv, d, A and the inequalities are given in 
the text. 
The solution of equation (17) is 14 



Hi = Z>^ + ^o) = exp[-( L: ^) 2 



He,ft + &) (26) 
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where D„(£ + £ ) is the parabolic cylinder function of order 17 and 
He, (£ + £ ) is the Hermite function of order 77. Only if a — * », 77 
becomes an integer, the Hennite function is reduced to a polynomial 
and Hi becomes the well-known solution of the parabolic lens-like 
medium extending to infinity. 3 

Outside of the guide, that is for r > R, the refractive index n is 
uniform, equation (14), and the solution of the wave equation (16) is 13 
the Hankel function of order v and argument kn r. That is 

H = Hl l) (kn r). (27) 

To match fields at the boundary r = R, the radial admittance 
H u /Ee inside and outside the guide must be identical. With the help 
of equations (15), (26) and (27), we obtain the characteristic equation 

. w nD&o) _ Hl l \knjl) (m 

k 2 D'M ~ H^XknoR) K Q) 

in which the derivatives are taken with respect to the arguments of the 
functions. 

We should have another boundary equation for the other side of the 
guide, r = R — a, but we are interested in guides with radius of curva- 
ture R small enough to push the field away from the center of the 
guide, and consequently the field at the interface r = R — a is negligibly 
small. 

To solve explicitly the boundary or characteristic equation (28) for 
/c 2 , we need asymptotic expansions of the functions involved. From the 
inequalities in equation (25), it can be deduced that 

I & I » 1 and I & I » b |. ( 29 ) 

The asymptotic expansion for D,(£o) is then 14 

DM = 6 »p (-| -in) + ^) S - »p (f) (30) 

where T(—y) is the gaussian function of argument (—»?)■ 

The asymptotic expansion for the Hankel function results from 
observing that as a consequence of equation (25) 

fcn ff B» 1, 

k : R»l, (31) 

^2^, 1 

k. 
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and 

(K - «)*p»i. 

Therefore we can replace the Hankel function by Watson's approxi- 
mation. 14 This approximation involves Bessel functions of order one- 
third and large arguments. Keeping the first term of their asymptotic 
expansions, the Hankel function results 

tf 1 ^ = aSS^' exp \& {kl - k M 

+ | exp[-^(fc» -«)*]}■ (32) 

Substituting equations (30) and (32) in equation (28) we obtain a 

simplified version of the characteristic equation 



V2t „..„,, _ (£_ 



l- r 7^)^exp^ + ^ 

1 + i exp [-|| {k\ - «)*] 

= *° (kl - fcV.)*w " (33) 

To solve this equation for k t we rewrite it as 

T(-v) = F( v ) (34) 

and notice that F(rj) is a large quantity. Therefore the gamma func- 
tion is also large and hence rj must be near a pole, which makes rj close 
to an integer p. Then we can replace the gamma function by the first 
term of the Laurent series (— 1) /p !(p — v)> and equation (34) becomes 

'->-£& (35) 

Substituting tj by the value given in equation (20) we derive the ex- 
plicit value of k, . This propagation constant is complex, k x = /3 + ia, 
and the real and imaginary parts are the phase and attenuation con- 
stants of the pth mode: 

/3 = Re k. 



=4-(^M + M i -i)+^]} 



(36) 
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a = Im k„ = 
where 



T (Rf l- d \ 3 ~\ 

exp L"3 vir) . i 



+ 2M - M 2 



dKRV2A 



(1 + M) 



(37) 



M = 



1 + 



f -\ p + l)id 



-h 



K = V2ttp\ 



(1-d) 2 

exp[f(l-^] 
[\/id(l - d)] 2p+1 



(38) 



(39) 



In equation (37), M affects the value of a mostly via the exponential 
and not via the fraction 

1 + 2M - M 2 



(1 + M) 

which for all practical purposes can be replaced by 1. Consequently 
the normalized loss per radian \/2ARa results 



L = V2A Ra m 



[ (Rfl- d\ 3 ~\ 

exp L" 3 \-ir) . 



dK 



(40) 



Now we turn to guides in which the refractive index is a function 
of y, Figs. 2a and 2b. 

Let us start with the ribbon-like structure of Fig. 2a and assume as 
in Ref . 6 that 



A, » A. 



(41) 



Provided that most of the electromagnetic field travels within the 
ribbon, the attenuation per radian is still given by equation (40), but 
the phase constant is a slight modification of equation (36). From Ref. 12 
is deduced 



A - ' - T 



*{q + 1) 
b 



2(1 - btfAiY' for field 



i+!t 



wb 

-2 



polarized along y, 

for field 
polarized along x, 



(42) 



where q + 1 indicates the number of maxima of electric field within 
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the guide along y and 

.1, - — £=• ( 43 ) 

Consider another guide, Fig. 2b, with rectangular cross-section and 
truncated parabolic index profile along both the x and y directions 



, =.[>- 4 +2 ^)' -Jgf] 



(44) 



Provided that most of the electromagnetic field is within the guide 
cross-section, the loss per radian is still given by equation (40), but 
the phase constant becomes 4 



i-li + i 1 



5 + 2 + 2 



-1 



vMSTM'^ni 



(45) 



M2 H wlkn 

where q + 1 is the number of maxima of the electric field along y and 



^2 = \/ 7=- (46) 

If 

•p = q = 

and 

a = b 

the guide has square cross-section and equations (40) and (45) yield 
a first approximation of the phase and attenuation constants in a 
curved SELFOC 9 guide. 
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